General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 




(NASA-*CR-I730 78) HOOELLING 

OF AOROBAL 

1583-35543^ 

ELBCTBOOEHAHICAL PROCESSES: 

□ HAGMETOSPHEBE 


TO MESOSPHERE Final Report 

(Acrcspace 


Corp* « El Sequndo, Calif.) 

26 p 

Uncial 

HC A03/MP AOI 

CSCL 04A G3/46 

36652 



A 


1 


• s 


✓ 


% 




t 




V 




t ^ ■ 


« 


r 




space Sciences Laboratory 
Report No. SSL-83 (7951)-4 


Final Re;iort 


Modelling of Auroral Electrodynamical Processes: 


Magnetosphere to Mesosphere • 


NASA Solar-Terrestrial Theory Contract NASW-3434 


Principal Investigator: 


Y. T. Chiu 

Co-investigators : 

D . J . Gorney 
A. M. Kishi 
A. L* Newman 
M. Schulz 

R. L. Walters cheid 
J. M. Cornwall 


S • S • Prasad 


Space Sciences Laboratory 
The Aerospace Corporation 
2350 E . ElSegundo Blvd 
El Segundo, CA 90245 


also at: Physics Department 

University of California 
Los Angeles, CA 90024 

Jet Propulsion Laboratory 
4800 Oak Groye Drive 
Pasadena, CA 91103 


Table of Contents 


Page 

Abstract ill 

X. Introduction ••• 1 

XX t Accomplishments of Program 3 

XXX . Discussion and Overview 15 

Appendix X: Publications and Presentations ...... 17 


PRECEDING PAGE BLANK NOT FILMED 


ill 


Abstract 


An overview of a three-year term of researches on auroral electrodynami- 
cal coupling between the magnetosphere and ionosphere-atmosphere is given o 
This final report succinctly summarizes the activities of the research team, 
which is comprised of scientists from The Aerospace Corporation, University of 
California at Los Angeles and Jet Propulsion Laboratory and is supported by 
the NASA Solar-Terrestrial Theory Program contract NASW-3434. It is shown 
that the interdisciplinary team has successfully focused its work upon 
building the basic elements in unifying local plasma theories of auroral 
magnetosphere-ionosphere coupling into a global one . The global goal of 
unification is being pursued in the new term of research. 
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INTRODUCTION 


This is the final report of the first three-year term of the Solar Ter- 
restrial Theory Grant entitled "Modelling of Auroral Electrodynamical Proces- 
ses; Magnetosphere to Mesosphere." 

In this report we shall attempt to give an overview of the direction and 
execution of researches performed under the contract * As the result of co- 
operative theoretical and observational researches , such as the present 
effort, scientific understanding of the auroral coupling between the magneto- 
sphere and atmosphere has grown by leaps and bounds in the last five or six 
years • A comprehensive assessment of this particular progress is documented 
in the magnetosphere-ionosphere working group report of the recent Solar- 
Terrestrial Physics Workshop at Coolfont, W. Va., and sponsored by NASA and 
NSF. The research program performed under the present contract is an integral 
part of the auroral work assessed in the above report; thus, the main 
viewpoint in this final report will be forward-looking, emphasizing what is 
ahead for auroral physics and how it fits into the global scheme of 
quantitative solar-terrestrial theory. 

We shall shiiw that the results of our researches under this contract are 
partial elements, which, together with results of auroral theoretical and 
observational researches of other groups, form a basis for unifying various 
theories of auroral magnetosphere-atmosphere coupling into a global scale 
theory • The unifying thread is the approach that we have used in our re- 
searches; global-scale plasma kinetic theory in which the global-scale mor- 
phology of magnetic and electric fields is the major driver of microscopic 
kinetic plasma phenomena which f orti the final link in the coupling . We have 
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successfully applied this point of view to formation of the evening discrete 

ft 

y arc, to the raorningside arc system and to polar cap arcs. With further obser- 
vational support from Dynamic Explorer and S3-3 results, we are now poised to 
attack the unifying problem of global arc formation theory. Indeed, as we are 
preparing this report, theoretical clarification of the relationship between 
polar cap arcs (including the recently discovered 6-auroral configuration) and 
oval arcs based on this point of view is at hand [Workshop on Polar Cap Pro- 
cesses, AFGL, Bedford, MA, 1983] . 

Aside from the forward-looking emphasis on unification of global auroral 
magnetosphere-atmosphere coupling, a second charac*>eristic of our researches 
under this contract is the emphasis on theory-data verifications and theoreti- 
cal guidance for accompanying observational programs . The close relationship 
of theoretical guidance and observational verification is particularly Weil- 
ls documented in S3-3 data analysis, from which an important portion of the 

j 

advances in auroral physics of the past 5-6 years has originated . An impor- 
tant achievement of our program is to carry forward this symbiotic relation- 
ship to understanding the kinetic plasma convection phenomena at equatorial 
' regions of auroral field lines as measured by the SCATHA satellite • Similar- 
ly, complex computation programs on auroral particle interactions with the 
atmosphere patiently built up under this contract will find their full verifi- 
cation in the theory of auroral emissions in the UV to X-ray range, which will 
be emphasized in future programs such as OPEN. The same can also be expected 
in the elucidation of dynamical atmospheric motions associated with auroral 
substorms . 

In, 

¥ . 
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II. ACCOMPLISHMENTS OF PROGRAM 


Since commencement of the present contract, in October 1980, researches 
undertaken in our program have progressed according to the plan envisioned in 
the initial proposal, A list of publications and presentations with acknow- 
ledged support, fully or partially, by the Sola t'-Terres trial Xheoi'y Program is 
given in Appendix I, 

In the initial proposal, the program plan was to Investigate a number of 
separate , but coherently related , topics by appropriate team members • At a 
mature stage of the investigation, these separate projects are to form the 
foundation from which a global auroral mqgnetosphere-ionosphere-atmosphere 
coupling model is to be constructed, i?br a summary of the physical basis for 
this staged development of the coupling model, the reader is referred to 
Section III of a review of the present status of auroral plasma research [Chiu 
et al,, 1983a] a This review also summarizes the cross-fertilization and 
cooperation with auroral experimenters, which characterise our work. Thus, a 
tangible accomplishment of our program is the guidance that our theoretical 
program provides to experimenters concerning the important physical issues in 
data analyses and plasma observations . These include developments in the 
observational study of auroral electric potential structure, ion conics, 
plasma injection into the auroral region and auroral return currents. 

In the past three years, several general areas pertaining to auroral 
magnetosphere-ionosphere-atmosphere coupling and forming separate elements of 
the global theory have been under investigation • The aecomplishments achieved 
in these areas are briefly reviewed below in terms of progress toward our 
proposed goal of developing a model of global magnetosphere-ionosphere- 
atmosphere coupling . 
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A» Electric Potential Structure in the Evening Sector 


A key element of auroral magnetosphere-ionosphere coupling is the 
understanding of discrete arc formation in the evening sector. The first 
efforts in explaining auroral E^ naturally concerned themselves with a local 
point of views a single auroral field line, with boundary conditions at the 
equator and ionosphere [Chiu and Schulz, J . Geophys . Res 83. 629, 1978] , 
This approach led to a fruitful understanding of how magnetic-mirror forces 
s^tt up E|j fields and roughly what strength these fields would have, but was 
not in principle capable of precise determination of auroral potential drops . 


It soon became clear that the ionosphere pJ.ayed an all-important 
role in coupling neighboring field lines [Chiu and Cornwall, 1980 j Lyons, J . 
Geophys . Res . . 86 , 17, 1980J , which thus forced an enlargement of the purely 
local point of view . An extreiuely valuable result of thit^ xater work was an 
understanding of what we will call the outer (or largest) scale of inverted- 
V's, that is, their latitudinal width. This width A is equal to (2p/Q)^/^, 
where is the auroral height-integrated Pedersen conductivity (as enhanced 
by electron precipitation) and Q •• n e^/mg where n is the density and Vg a 
typical velocity of plasma-sheet electrons which will be accelerated by the 
auroral potential drop . It turns out that A * 50 - 100 km, in excellent 
agreement with experimental data . 


This outer scale links the electric fields along S with those across 
and thus the global morphology of electric and magnetic fields [Lyons, J . 
Geophys, Res. , 86 , 17, 1980], according to a formula (Chiu et al., 1981] which 
we quote only in a highly-simplified form: 


A(J) “ - / dx e ^ ^ E,(x) 
0 ^ 


( 1 ) 
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where x is N-*S distance in the ionosphere, and is the equatorial perpen~ 
dicular field mapped to the ionosphere as if there was no parallel potential 
drop. The left hand side, A(J>, is the (positive) potential drop between iono- 
sphere and equator along the center (x-0) of the inverted-V, and in this 
simple form we take it that Ej^(x) " Ej^(-x) , that i<s, there is a reversal of 
the perpendicular convection field about the central auroral field line. 

Thus, at one stroke, we are able to bridge the gap between local 
one- and two-dimensional theories to global two-dimensional theories in which 
the global E^^ structure is important . 

Observational verification of the various aspects of (1) indicates 
to us that our basic kinetic approach to evening auroral magnetosphere- 
ionosphere coupling is on the right track. Meanwhile, other key elements 
toward a comprehensive model of the evening electric potential structure have 
been investigated. These are: (1) elucidation of the relationship between 
observed adiabatic plasma motion and plasma turbulence effects, and (2) solu- 
tion of Poisson’s equation for the two-dimensional structure of which is 
considerably more complicated than the qtiantity A$, the ionosphere-to- 
magnetosphere potential drop • Researches on both areas are progressing 
according to schedule and interim results have been, or are being, pub- 
lished. Cornwall and Chiu [1982J, in particular, have shown that plasma 
turbulence can have very important time-dependent effects on the global struc- 
ture of the coupling. In this work electrostatic ion cyclotron waves, driven 
by an electron beam, heat ions in T^^ thus reducing the magnetic-mirror forces 
which set up the Ep which drives the electron beam. The result can be a 
steady state accommodating of both adiabatic and turbulent effects, as 
observed, but it can also be a relaxation oscillator whose longer time period 
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is •» 10 sec, identified with an ion-transit time scale 


Preliminary results of a detailed two-dimensional auroral potential 
structure model, involving complex equilibria of many species of auroral 
plasma and static electrostatic field#, have also been communicated to the 
space plasma physics community [e.g. presentation by Newman et al., 1982 J . A 
comprehensive paper on this topic is currently in preparation [Newman et al., 
1983] • A significant theoretical result of this complex model is the proposed 
description of the auroral return current as driven by a downward parallel 
potential drop of 100 eV thus driving ionospheric electrons upward to 
balance the precipitaticn of magnetdspheric electrons in the auroral arc. 
This proposal, originally advanced in Chiu, Newman and Cornwall (1981), has 
turned out to be prophetic in that evidence of upward eleesron beams of 
similar energy has recently been discovered la Dy nami c Explorer results [Burch 
et al., Geophys . Res . Lett . , 10 , 753, 1983] and in S3-3 data. The S3-3 data 
further Indicates [Gorney, Croley and Chiu, 1983] that this parallel electric 
field plays an Important role in trapping ionospheric ions, thus allowing them 
to be heated to high energies to form conics • This may be the initial step in 
arriving at the solution to the long standing puzzle of conic formation. 

B . Morningslde and Dayslde Auroras 

To construct a global model of auroral magnetosphere-ionosphere- 
atmosphere coupling, an understanding of aurora formation mechanisms in the 
morning and dayslde sectors are just as important as that of the evening 
sector, even though the current theoretical effort in the scientific community 
is concentrated on the parallel electric field of the evening sector. Al- 
though observations of auroral plasmas have long indicated absence of electric 
acceleration in the morningside auroras, a systematic attack on their forma- 
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tion raechatiiam has not been attempted# By invoking recently-analyzed results 
of plasma observations on S3-3 and SCATHA satellites ^ we have constructed a 
theory of morningside arc systems based on the modulation of diffuse auroral 
electron precipitation by the mirror instability driven by auroral ions [Chiu 
et al#, 1983b]# An especially interesting feature of th'^s theory, in agree- 
ment with SCATHA plasma observations, is that the ionospheric ions (0***) accel- 
erated upward from the evening discrete arc and convected to the diffuse 
aurora region is a key element in the stabilization of the mirror instability 
in the evening sector; thus, again demonstrating the necessity of considering 
auroral magnetosphere-ionosphere coupling in its global setting • 

C # Polar Cap Auroras 

Correlation studies of polar cap arcs and low-altitude polar cap 
convection patterns (Ej^) have indicated a relationship between these integral 
elements of auroral magnetosphere-ionosphere coupling to the pattern of north- 
ward interplanetary magnetic fields. Recently, the importance of such corre- 
lative studies has been brought into focus by the new discovery fo the 9-auro- 
ra [Frank et al#, Geophys # Rev # Lett 1001, 1982] for which the low- 
altitude electron signatures were noted to be that of the inverted-V. 

From the perspective of global aurora formation theory, the unique- 
ness and, at the same time, similarity of the polar cap aurora with other 
auroras pose a major question in auroral physics: Are we witnessing a new 
type of aurora formatitsn process entirely distinct from the evening discrete 
arc, with signature, and from the morningside arc system, without Ej sig- 
nature [Chiu et al#y 1983a]? 
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Chiu and Gorn^sy U983J dfesjons Crated that polar cap arcs are Identi- 
cal to oval discrete arcs in electric field and plasma signatures . Evidently , 
Intrusions or presence of hot plasma from the plasmasheet and/or magnetosheaCh 
produces discrete arcs of whatever configuration if Che proper signature 
indicating upward field-aligned current is also present. In short, arc forma- 
tion of whatever configuration takes place if (I) is satisfied. Thus, the 
combined consequence of (1) and thfr plasma Intrusion hypothesis is that arcs 
formed by parallel potential drops (oval and polar cap discrete arcs) of 
whatever configuration have their origin in the concentrated of the boun- 
dary layer, which, in whatever configuration, is also the source of the hot 
plasma being precipitated by the parallel potential drop. Moreover, this 
conclusion, though illustrated with the kinetic theory of arc formation above, 
has been demonstrated with plasma ani electric field data and is not limited 
hrt any one particular theory • 

The theme of global unification of arc formation theory is once 
again demonstrated. It now remains for our researches to advance the hypothe- 
sis by examining how auroral plasmas come to be present at auroral field lines 
and with what kinetic signatures . 

D . Auroral Plasma Formation 

A 

In the past few years, researches in magnetosphere-ionosphere coup- 
ling have focused considerable attention upon how auroral plasmas precipitate 
to form the optical signature. However, underlying the precipitation mecha- 
nisms, a major question of kinetic aurora theory needs to be addressed in 
detail; From where and how does the auroral plasma acquire its kinetic char- 
acteristics such as temperature and pitch angle distribution? In short, how 
is the auroral plasma formed? 
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By using Che propcrCies of the exact Green’s function solution of 
the collisionleoe time-dependent Boltzmann equation, Chiu and Kisht (1983) 
have been able to construct a time-dependent model of auroral plasma formation 
by global time-dependent convection electric fields in dipolar magnetic geo'Sm'^ 
try. The significant result of this development is that two-dimensional 
(V||,Vj^) distribution functions can now be mapped backward and forward in 
space-time in an economical invariant-preserving manner, Since diS/tribution 
function information in Che equatorial auroral region is now made available in 
recent satellite particle experiments such as those onboard SCATHA and GEOS, 
we expect that a natural consequence of our development will be the ability to 
interpret the satellite distribution function data in terms of sources and 
fields; whereas, previously, one is restricted to segmented tests of the 
convection model such as energy-time dispersion in Isolation of Important 
parameter^? such as pitch-angle distribution. 

Because of the large variety of theoretical and observational appli- 
cation of our model, we have restricted Che presentation in the above initial 
work to the time-dependent effects upon the convection of electrons. Even so, 
we have demonstrated that many of the interesting features of electron 
distribution functions seen for the first time recently can be schematically 
understood in terms of a very simple time-dependent convection model , We have 
shown examples that bear on the formation of square-box gradients and field- 
aligned dlstribu?;ions, on injection fronts and on inverse dispersion signa- 
tures. The principles of detailed data verification of theory, using SCATHA 
data, will be published in cooperation with experimenters . 
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E. Electrodynamlc Coupling with the Thermosphere 


The modeling of the dynamical response of the high latitude neutral 
thermosphere to sudden deposition of auroral energy, using an idealisied high- 
resolution spectral model, is advancing according to the plan stated in the 
initial proposal. The model [Walterscheld and Boucher, 1982 J treats distur- 
bances of the thermosphere due either to the deposition of energy by a momen- 
tum source (such as ion drag) or a heat source (such as Joule dissipation or 
particle precipitation) . The model describes two-dimensional dissipationless 
motion on a basic state which is isothermal and at rest. The model differs 
from steady-state models (such as employed by investigators at Utah State 
University) ir that it is time dependent, and differs from previous idealized 
time-dependcrnt models in that it includes the effects of the earth’s rota- 
tion, Rotation permits the existence of inertial oscillations and geostroph- 
ically balanced motion. The model is appropriate for studying the transient 
response of the lower thermosphere to typical auroral and sub-auroral scale 
disturbances . These scales may be considerably smaller than those of present 
global models such as at NCAR and GSFC , Initial calculations were performed 
for selected vertical wave numbers in order to study the dependence of the 
dynamical response on the vertical scale of the motion, ' Model results indi- 
cate that: (1) for an ion drag source a strong residual geostrophic flow is 
favored by forcing in a deep layer, whereas a strong residual inertial oscil- 
lation is favored by forcing in a shallow layer, (2) for a Joule heating 
source a strong residual geostrophic flow and inertial oscillation are both 
favored by forcing in a shallow layer, and (3) the ion drag source may excite 
gravity waves through the action of the Coriolis force in deforming pressure 
surfaces across the initial flow, the effect being more pronounced for forcing 
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in a shallow layer • We have extended our model to consider contributions from 
a continuous spectrum of vertical scales of motion excited by an Initial 
disturbance of a prescribed width and depth. The results were in agreement 
with the findings of the earlier analysis • A somewhat surprising result was 
that even for comparatively deep forcing a non-negligible residual flow 
evolved for the heating case. This was due to the fact that, a given pertur- 
bation in the height of a pressure surface, shallow oscillations will generate 
larger velocity perturbations than deep oscillations, thus the effect of the 
shorter scales is enhanced relative to their prevalence. 

Our initial modeling effort 'involved time-dependent calculations of 
the transient response of the atmosphere to a sudden injection of energy and 
momentum. A simple linear model Indicated that a rapid deformation of 
pressure surfaces across an initial eastward flow can occur (unless the 
initial disturbance is very shallow), bringing the flow into geostrophic 
balance with little reduction in amplitude. However, if the initial flow is 
very shallow, adjustment takes place slowly, and much of the energy of the 
initial flow is transferred to inertia-gravity wave. A paper on the subject 
has been revised and will be published in the Journal of the Atmospheric 
Sciences « In addition, an invited review paper on the subject of geostrophic 
adjustment was presented at the Spring Meeting of the Merican Geophysical 
Union . 

In collaboration with Prof . Karl Taylor of the University of 
Florida, who is spending the summer working in the laboratory, we have begun 
development of a sophisticated high-resolution two-dimensional model of the 
atmospheric response to a geomagnetic substorm. Our goals are to study in 
detail the dynamical adjustment to time-dependent heat and momentum sources 
when dissipative and nonlinear effects are taken into account, and to develop 
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nonreflecting boundary conditions for use in a high-resolution three- 
dimensional model* (Since a high-resolution three-dimensional model would be 
prohibitively expensive for a large domain of integration, boundary conditions 
must be implemented which will prevent fast moving waves generated in the 
source region from being reflected back into the auroral region during the 
simulation.) We have formulated the model equations, worked out a numerical 
scheme, substantially completed coding of a model with homogeneous boundary 
conditions , and have made some initial runs , 

F , Auroral Electron Interaction with the Atmosphere 

The stated problem was studied with respect to: (a) the high energy 
(0.5-5 keV) electron precipitation occurring in the discrete auroral arcs of 
the evening sector, and (b) with respect to the much softer electron (a few 
teais of eV - iOO ev) precipitation taking place in the cusp region • 

Interaction of auroral electrons with the atmosphere produces 
ionization and excitation of atmospheric species and a significant 
backscattering of both primary and secondary electrons in the upward 

direction. Backscattered electrons with E < Eb trapped between the lossy 
ionosphere below and the electric field above, and are reflected back and 
forth a few times before precipitation into the atmosphere . This important 
effect was Included in our studies of auroral electron interaction with the 
atmosphere . We used numerical solutions of the multiangle equation of 
transfer at all energies. Most previous studies of auroral electron- 
atmospbiare interaction employed a hybrid approach which Uses Foi\ker*^Planck 
techniques in the range E > 500 eV and transport equation at lower >inergies . 
This approach obviously suffers from the problem of matching the solution at 
the boundary of the two energy regions . We avoided this problem by solving 
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the transport equation at all energies • Another serious shortcoming of all 
previous studies of this problem has been that they considered only two 
streams, vix., an upward beam and a downward beam. In contrast, we considered 
sixteen pitch angles, eight corresponding to the upward directions and an 
equal number in the downward direction. Due to this significant improvement 
we were able to obtain the pitch angle distribution of the auroral electron 
throughout the atmosphere. As a result we found that the auroral plasma, 
including those observed at the satellite and rocket altitudes, has three 
components: (a) a low energy isotropic component consisting of secondary 
electrons of atmospheric origin, (b) mostly downward flux which included the 
auroral primary beam, and (c) a field perpendicular component with energy 
around the beam energy . The identification of a field perpendicular component 
constitutes an important outcome of the present study, because the existence 
of this component can lead to interesting kinetic effects , 

Interaction of the softer electrons of the cusp region with the 
atmosphere presents a more difficult problem with respect to the solution of 
the transfer equation. The elastic scattering of high energy auroral elec- 
trons can be satisfactorily described by the Rutherford cross section . Such 
is not the case with the low energy cusp electrons . The differential elastic 
scattering cross sections of low enei'gy electrons (E < 20D eV) tend to have a 
peak in the backward direction . This becomes quite important at or below 100 
eV. The transport code used for the auroral study was therefore modified to 
take into account this feature in the differential elastic scattering cross 
sections of low energy electrons. For this .cpose, the measured differential 
elastic scattering cross sections of 0, O 2 and N 2 were fitted by a sum of two 
Rutherford type expressions . The original code used a numerical method to 
evaluate certain integral over pitch angles . The modified code performed the 
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same integral * 3 nalytically • This improvement jts quite useful in the study of 
the cusx^ region, because the differential elastic scattering cross section of 
low energy electrons involved two Rutherford type cross sections instead of 
one* Another special feature of the cusp region is the role of photoelec** 
trons* In our studies of the cusp region, therefore, we examined the fluxes 
of both photoelectrons and precipitating soft electrons. Depending upon the 
time of the day (i.e?, solar zenith angle) the fluxes of the photoelectrons 
may completely swamp the fluxes due to the precipitating electrons in the low 
energy (E < 20-30 eV) region. The photoelectron population, however, 
diminishes dramatically above 50-60 eV. Observations in this high energy 
region, therefore, are representative of the precipitating electrons 
irrespective of the solar zenith angle. At twilight (solar zenith angle > 
100®), the contribution of the photoelectrons is, of course, negligible. 

Currently, at the completion of the above studies, we are directing 
our attention towards application of the atmospheric interaction code to 
calculate the emission of auroral X-rays and UV electromagnetic radiation. 
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III. DISCUSSION AND OVERVIEW 


Research activities reviewed in the previous section and publications 
listed in Appevldix I summarize the execution of our research plan, as envi- 
sioned in the original proposal in 197? ^ In sum, the research team, composed 
of both senior and junior scientists of several distinct solar-terrestrial 
disciplines, channeled their activities toward a single goal: understanding 
of auroral magnetosphere-ionosphere-atmosphere coupling . Included in our 
research activities are frequent visits and cooperative work with scientists 
from other institutions • 

Since we are at the base-building phase of this long-term goal, our 
initial researches seem to cover a spectrum of auroral phenomena whose 
relationships tn each other have yet to be fully explored . Such a phase of 
research activity is unavoidable if our research plan is to be realized . This 
is especially true since complex computational programs, such as modelling of 
auroral potential structures, auroral atmospheric dynamics and auroral elec- 
tron interactions with the atmosphere, must be individually constructed to 
form the basis for our eventual goal of unification. 

Despite the distraction of such necessary preparatory work, the physical 
content of our basic theme of seeking the origins of auroral magnetosphere- 
ionosphere-atmosphere coupling in global scale kinetic plasma theory is re- 
vealed in our emphasis upon observational verification of our theoretical 
results. To make this point clear, it is perhaps important to carefully 
distinguish our approach to auroral physics ( global scale kinetic plasma 
theory) from simulations and other local scale considerations of microphysics 
in which global scale effects are set up by artificial boundary conditions of 
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the simulation "box” • Our cooperative examination of auroral observations 
with experimentalists has demonstrated, again and again, that the two ap- 
proaches are not the same (c£. Chiu et al., 1983a). Global scale auroral 
plasmas, in its cumulative interactions with global scale inhomogeneous mag- 
netic and electric fields , arrives at the auroral zone with its own boundary 
conditions which In turn are influenced by the microphysics . Local microsco- 
pic interactions, amenable to simulation studies, are but one of the many 
factors that influence the behavior of the plasma in the auroral zone * Equa- 
tion (1) is a very good example of this interaction between "boundary condi- 
tions" and kinetic plasma characteristics • The parallel electric field has 
its origins in the perpendicular electric field which in turn is dictated by 
the global configurations of the magnetosphere . Microscopic plasma interac- 
tici'fis merely deals with the distribution of the parallel electric field . 
Thus, to seek the origin of the discrete arc is to seek the origin of the 
perpendicular electric configuration favorable to upward field-aligned cur- 
rent, as indicated in Eq . (1). Since simulation studies are (at present) 
incapable of dealing with these global issues, a research strategy based on 
the Boltzmann-Maxwell equations must be used to connect the spectrum of scales 
involved in auroral phenomena. The example of aisroral potential structure, 
cited above to illustrate this principle, is not unique. We find, for 
example, the same operating principle in SCATHA plasma data which show that 
global electric and magnetic field structures dictate the plasma kim^tic 
characteristics observed in the auroral equatorial region (Chiu and Kishi, 
1983). 

The unification of global auroral structure and auroral microphysics 
is precisely the theme of the next phase of our research plan • 
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